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Abstract 

Convection in horizontal layers of binary fluids heated from below and in 
particular the influence of the Soret effect on the bifurcation properties of 
extended stationary and traveling patterns that occur for negative Soret cou- 
pling is investigated theoretically. The fixed points corresponding to these 
two convection structures are determined for realistic boundary conditions 
with a many mode Galerkin scheme for temperature and concentration and 
an accurate one mode truncation of the velocity field. This solution proce- 
dure yields the stable and unstable solutions for all stationary and traveling 
patterns so that complete phase diagrams for the different convection types 
in typical binary liquid mixtures can easily be computed. Also the transition 
from weakly to strongly nonlinear states can be analyzed in detail. An inves- 
tigation of the concentration current and of the relevance of its constituents 
shows the way for a simplification of the mode representation of temperature 
and concentration field as well as for an analytically manageable few mode 
description. 

PACS number(s): 47.20.-k, 47.10.+g, 51.10.+y 
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I. INTRODUCTION 

Convection in binary miscible fluids like, e. g., ethanol-water or ^He-^He is a well estab- 
lished and accepted system for studying instabilities, bifurcations, complex spatiotemporal 
behaviour, and turbulence. This is on the one hand due to its sufficiently simple experimen- 
tal realization under well controllable conditions. On the other hand, a great advantage for 
the theoretical analysis is the solid knowledge of the governing field equations. So, recently a 
lot of research activities ||T]-0 have been directed towards investigating the enormous vari- 
ety of pattern forming behaviour in this system. The richness of spatio temporal phenomena 
in binary fluid mixtures stems from a feed-back loop between the fields of velocity, concen- 
tration, and temperature. Let us start with the velocity field: The convective flow is driven 
by the buoyancy force field which itself is determined by variations of the temperature and 
of the concentration field. The latter are on the one hand generated via the thermodiffusive 
Soret effect by temperature gradients and on the other hand they are reduced by concentra- 
tion diffusion and by mixing due to the convective flow. Since these changes influence the 
buoyancy which drives the flow the feed back loop is closed. 

In this article we concentrate on the Soret coupling and its influence on spatially extended 
convection states of straight parallel rolls that occur either as a horizontally traveling wave 
(TW) or in the form of stationary "overturning" convection (SOC) rolls. Among others, we 
elucidate the Soret induced changes in the combined SOC-TW bifurcation topology which 
offers in both types of convection the possibility of sub- and supercritically bifurcating 
branches depending on the strength of the Soret coupling. Both solution branches develop 
with increasing Soret coupling saddle node bifurcations which give rise to stability changes. 
Finally, there exists a merging point of the SOC and the TW branches for moderate negative 
Soret couplings. This competition of stationary and traveling states can only be observed 
for negative Soret coupling, where temperature gradients induce adverse concentration gra- 
dients that stabilize the unstable thermal layering. For positive Soret couplings, there is no 
oscillatory instability of the basic state. The two interesting cases for negative coupling are 
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now weak and strong Soret effect. For the latter one a bistability of slow and fast TWs was 



recently reported [|T^ which coexist with the likewise stable basic state of heat conduction. 



A detailed phase diagram was discussed |T^ describing the Soret dependence of the saddle 
nodes. For weak Soret couplings, however, a detailed study of the bifurcation toplogy was 
missing. Data from direct numerical simulations are sparse since in the vicinity of saddle 
nodes and bifurcation points the intrinsic time scale of the system is arbitrarily long. Thus, 
phase diagrams to elucidate the whole bifurcation topology were incomplete. 

We have determined the SOC and TW fixed points of the system by a many mode 
Galerkin scheme whose convergence properties do not depend on the time scale of the sys- 
tem. In particular, this allowed us to study all unstable branches on which, especially, the 
transition from weakly to strongly nonlinear convection takes place. Furthermore, a detailed 
explanation of the concentration distribution and its relation to convection is given. 

Our article is organized as follows: The second section presents the system, the fields 
needed for its description, their governing equations with the explanation of the relevant fluid 
and control parameters. Finally, it presents a short survey on the typical bifurcation scenario 
in the convection of binary liquid mixtures with negative Soret coupling. The third section 
shows the field truncations, the method of solving the system and the solutions of our many 
mode Galerkin scheme basing on a reasonable approximation in the velocity field. In the 
fourth section we discuss the influence of the Soret coupling on the bifurcation topology 
by means of exemplary bifurcation diagrams realized in experimentally feasible mixtures 
and a detailed phase diagram. Furthermore, evidence for an instability of a TW towards a 
modulated TW (MTW) is given. The fluid parameter range for its occurence is elucidated 
for ethanol-water as well as for ^He-^He-mixtures. Finally, we extract the relation between 
concentration distribution and convective structure and we investigate the importance of 
the Soret effect at the boundaries and its neglibility in the bulk. 
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II. SYSTEM 

A layer of a binary fluid mixture with a mean temperature T and a mean concentration 
C is confined between two perfectly heat conducting and impervious plates separated by a 
distance d and exposed to a vertical, homogeneous gravitational acceleration g. The lower 
(upper) plate is kept at a fixed temperature T + AT/2 (T — AT/2). 

The fluid parameters are p (density of the fluid), = (thermal expansion co- 

efficient), P = —^§^ (solutal expansion coefficient), u (kinematic viscosity), k (thermal 
diffusivity) , kx (thermodiffusion coefficient), and D (solutal diff usivity) . 



A. Scaling and balance equations 



We scale lengths by the height d of the layer, times by the vertical diffusion time d'^/n 
of the heat and accordingly velocities by n/d. The deviation T of the temperature from its 
mean T is reduced by AT, that of the concentration field by ^AT, and the pressure p by 
Then, the balance equations for mass, momentum, heat, and concentration [Q,^ read 



in Oberbeck-Boussinesq approximation |18 



dtu 



V - u 

-(u- V)u- V 



p + 



d^ 



9 z 



+aV^u + Ra{T + C) e, 
dtT = -V -Q = -(u ■ V)T + V^T 
dtC = -V -3 = -(u • V)C + LV^ {C - ipT) . 



(2.1a) 

(2.1b) 
(2.1c) 
(2.1d) 



The Dufour effect describing currents of heat driven by concentration gradients is discarded 
in ( p.lcp since it is relevant only in binary gas mixtures [0,|1^] or in liquids near the critical 
point [|g. 

The dimensionless fluid parameters are the Prandtl number a = u/k, the Lewis number 
L = D / K, and the separation ratio ip = — f The latter characterizes the strength of the 
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Soret effect. The Rayleigh number R = serves as control parameter measuring the 

thermal stress. 

B. Bounday conditions 

We use experimentally realized boundary conditions for the top and bottom plates at 
z — ±1/2 which are no slip for the velocity field, 

u(x,y,z = ±l/2;i) = , 

perfectly heat conducting for the temperature field, 

T{x,y,z = ±1/2- t) = Tl/'2 , 

and impermeable for the concentration field, i. e. 

e,-J = -Ld,{C-il:T){x,y,z = ±l/2;t) = . (2.2) 

We restrict ourselves to the description of extended roll like patterns that are homogeneous 
in one lateral direction say y. So, we investigate two dimensional states of a certain lateral 
periodicity length X — 2Tr/k. In most cases we take A; = tt, i. e., A twice the thickness of 
the fluid layer, which is close to the critical wavelengths for the negative Soret couplings 
investigated here. Furthermore, the stable nonlinear TW and SOC states that are observed 
in experiments have typically a wave number k — tt. 

C. Conductive state 

In the motionless basic state, a vertically linear temperature profile, Tcond = —z, is 
observed due to the different top and bottom temperature. This leads via the Soret effect 
and the no flux condition for J to a likewise linear concentration profile, Ccond = —ipz. Both 
together yield the hydrostatic pressure 

1 . . 2 fd^ \ 
Pcond = Po - + V')^ - I -^9 1 ^ 

in the quiescent state. 
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D. Control and order parameters 

The dimensionless temperature difference between tlie two plates, namely tlie Rayleigfi 
number R, is used as contol parameter. Mostly we scale it by the value of the onset of 
convection in a pure fluid: 

- A - ^ 
^ ~R0- 1707.762 ■ 

The convective states of the system are characterized by four order parameters: 

(i) The maximum Wmax of the vertical velocity field. 

(ii) The Nusselt number = (Q ■ e^)^ giving the lateral average of the vertical heat 
current through the system. In the basic state of heat conduction its value is 1 and larger 
than 1 in all convective states. 

(iii) The variance M = a/ {C'^)x,z/ {Ccond)^,z concentration field being a measure 
for the mixing in the system. The better the fluid is mixed the more the concentration is 
globally equilibrated to its mean value so that M vanishes in optimally mixing, strongly 
convecting states. 

(iv) The frequency ci; of a traveling wave. Thus, extended TWs with a wave number k 
have a phase velocity v = uo/k. They are stationary states in a reference frame comoving 
with V relative to the laboratory system. 

E. Typical bifurcation scenario 

For fluid parameters typically realized in mixtures of water and about 10 wt.% ethanol, 
an oscillatory, subcritical onset of convection is observed. It is connected by an unstable 
TW branch with a saddle node bifurcation giving rise to stable, strongly nonlinear TW 
states. At a certain Rayleigh number, the phase velocity of these waves vanishes and the 
SOC branch of stable stationary states can be observed. Along the TW bifurcation branch 
which is shown in detail in Fig. the concentration changes its structure from lateral ho- 
mogeneity and vertically linear layering in the basic state over plateau-like distributions in 
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fast TWs to boundary layer dominated, slowly traveling waves (see the discussion and the 
figures in Sec. |VB|). The contrast between two adjacent TW rolls is strongly related to the 
phase velocity of the TW and it vanishes with this velocity. So, SOCs do not show such a 
concentration contrast. In SOCs adjacent rolls are mirror images of each other and they are 
separated from another and from the top and bottom plate only by thin boundary layers. 
The latter are a typical phenomenon for convection of weakly diffusing scalars. 

III. COMPUTATION OF EXTENDED STATES 
A. Modelling the velocity field 

In liquid binary mixtures like ethanol- water momentum diffuses approximately ten times 
faster than heat. This means that the Prandtl number a is of 0(10) so that the velocity field 
may be adiabatically eliminated. Then, the momentum balance ( |2.1b| ), say in vertical direc- 
tion, reduces to the balance of the diffusive term aV^w and the buoyant term Ra (T + C) 
with the latter containing no derivatives. Thus, in a stationary flow, either in the laboratory 
frame or in a comoving one, the amplitudes of higher lateral Fourier modes Wn of the vertical 
velocity fleld w = u ■ scale at least like so that they decrease rapidly and even faster 
than those of the temperature fleld. That is the reason why higher modes than the critical 
flrst lateral Fourier mode are not necessary for a good description of the velocity fleld. One 
can expect this to hold for all a > 1. 

The next question deals with the role of the lateral mean of the velocity fleld, i. e., 
of its zeroth lateral Fourier mode. Continuity implies that only the lateral velocity fleld 
M = u ■ Gj. can contain such a mean flow. In order to determine the relevance of a z- 
dependent mean flow we compare its maximum with the two other velocities in the system: 
The maximal vertical flow velocity Wmax and the phase velocity v of TWs. Just at the onset 
of convection, the mean flow may be estimated ||21]] to scale with w^ax ^"^^ be very small 
in comparison with v. Furthermore, mean flow and phase velocity have the same sign. In 
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strongly nonlinear states, the ratio of flow velocity and mean flow is nearly 10^ |T^, but 
has changed sign. Thus, the mean flow has a non-monotonous dependence on the phase 
and flow velocity whereas all other properties of TWs like mixing and heat transport vary 
monotonically. Hence, the mean flow cannot contribute systematically to these properties 
characterizing TWs sufficiently. This is the reason, together with the smallness of the mean 
flow, for ignoring it altogether. 

The z-dependence of the critical velocity field is described in an adequate manner (see. 



e. g., [^,^) by the first even Chandrasekhar function Ci{z). Then, the velocity field 
of straight rolls with axes oriented in y-direction that are propagating with phase velocity 
V in x-direction is described by 

/ \ 

— i sin k{x—vt)C[{z) 



u{x, z; t) 



Ci(0) 



v 



COS k{x—vt) Ci{z) 



(3.1) 



Herein, the phase is chosen so that the maximal vertical flow occurs at t = and x = 0. 

Fig. checks in the first row the applicability of the ansatz ( p.l| ) by plotting the contribu- 
tion of modes in the vertical velocity w which are not represented by ( |3.1| ) for two separation 
ratios ip. The convective amplitude t^max was chosen as abscissa to quantify the nonlinearity 
of the states. As a measure for the strength of higher contributions to the vertical velocity 
the error 



{^MAc)x,z - {wWmAc)x,. 
{^MAc)^,z + {wWmAc)x,^ 



1/2 



(3.2) 



was computed where WMAci^y^) denotes the velocity field calculated from the full field 
equations by means of a finite difference MAC scheme |T6|. w{x, z) is the one mode approx- 
imation ( |3.1| ) with the two velocities v and Wmax chosen such that the numerically obtained 
velocity WMAci^^z) is fitted best or, equivalently, the error (|3.2|) is minimized. For both 
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Soret coupling strengths shown in Fig. |^ the error of the ansatz is smaller than 4% but 
increases as expected with the convective amplitude. 

The bottom row of Fig. ^ shows the maximum of the mean flow {u{x,z))x- For both 
separation ratios the mean flow is smaller than the convective amplitude by at least three 
orders of magnitude. 

All in all, Fig. ^justifles the approximations implied by the ansatz for the veloc- 
ity fleld. An important consequence of the flxed spatial structure of the velocity fleld is 
that all nonlinear terms in the balance equations contain now the same amplitude, namely 
Wmax; slucc all nonlinearitics are convective ones. This is an enormous simpliflcation for the 
theoretical analysis as we will see below. 



B. Galerkin expansion for temperature and concentration 



The temperature fleld T{x, z; t) is appropriately described by 

oo oo 

T{x, z;t) = + [TiTi^y^^"'"^'" + c.c] ^2 sin2m7r;z 

n=0 m=l 
oo oo 

+ I] I] [r|r+V(t)e-*('"+^)'=" + C.C.] y2cos(2m + l\ 

n=0 m=0 

This representation incorporates the mirror glide symmetry 



TTZ 



T 



C 



> (x, z; t) 



T 



C 



of TW and SOC states 13,16 



(3.3) 



(3.4) 



The representation of the concentration fleld is more subtle because of its boundary 



condition ( p.2|) coupled to the temperature fleld. The straight forward solution is the intro- 
duction of the combined fleld 



C(x..;t) = -C(...;<)-r(...;<) 



(3.5) 
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obeying the equation 

dtC = -{u-V)C + LV\ - V'T (3.6) 
and the boundary condition 

d,C{x,z = ±l/2;t) =0 . (3.7) 



Note that the ( field (|3.5| ) differs by a scahng factor l/ifj from the field that has mostly been 
used so far, see, e. g., i25| , |2TlJl9(| . An adequate trigonometric expansion is 



oo oo 

({x, z,t) = J2Yl [C2r'(^)e-'^"'^ + C.C.] v^sin(2m + l)7rz 

n=0 m=0 

oo oo 

+ E E [C2r+i(i)e-^('"+'^'^ + c.c] y2cos2m7r^ (3.8) 

n=0 m=0 

which also takes the mirror glide symmetry ( p.4| ) into account. 

The introduction of the combined field C was also motivated by the wish to fulfill the 
boundary condition (|2.2|) or (|3.7|) for the concentration field exactly. However, the formula- 
tion of the concentration balance in terms of the ( field causes a severe theoretical drawback 
for small Lewis numbers L and separation ratios ^p of order 0(1), that are typically realized 
in liquid mixtures: ( and T have the same order of magnitude according to (|3.5|) . But in 
the balance equation (|3.6|) for ( the diffsuive term enters with weight L = 0(0.01) and the 
temperature with 0(1). This means that for an appropriate solution of (|3.6|) for a partic- 
ular (^-mode temperature modes are necessary which are 1/L times, i. e. about 100 times, 
smaller than the ^-mode under consideration. Despite the fact that a relevant contribution 
of higher temperature modes was not observed neither in experiments nor in simulations 
they are necessary in ( ^.6[ ) if the ^-field is introduced and small Lewis numbers are con- 
sidered. That is mainly the reason why earlier Galerkin approximations using the ( field 
and only few temperature modes did not succeed in describing nonlinear TW convection in 
binary liquid mixtures [plip3ip5H 27[] . In binary gas mixtures with typically L = 0(1) this 



problem does not occur [0 
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C. Solution of the system of mode equations 

After projecting the balance equations for T and ( onto the bases used in the mode 
expansions ( p.3|) for T and (|3.8|) for ( one gets ordinary differential equations for the mode 
amplitudes of the temperature and the ( field 

x(t) = [r^'f it) , Tir+v it) , cir ' it) , cir+i it)] 

The column vector X is written here as the transpose of a row vector. 

In the case of SOC the amplitudes are constant and have to be chosen real in order to 
be compatible with the velocity ansatz (|3.1|) for v = 0. The flow amplitude w^a.^ appears 
via (|3.1| ) in the T and ( field equations via the convective nonlinearity and in addition as an 
inhomogeneous contribution wdzTcond = —w from the conductive part, Tcond = —z, of the 
temperature field ( p.3|) . Thus, the mode equations for SOC states take the form 

) ■ X — Bsoc( ^max ) • (3.9) 

Here, the matrix AlgQ(^(^max) of mode coupling coefficients contains if max from the con- 



vective nonlinearity. The momentum balance in (|2.1b| ) provides the relation between the 
Rayleigh number R and the velocity amplitude Wmax 

R= ^ (3.10) 

a - Xsoc(ti^max) 

containing the solution Xsoc = -M^}. ■ Bsoc of ( p.9|) depending on Wmax- The vector 
a contains all projection coefficients. The pressure gradient in the momentum balance of 
( p.lbp may be eliminated by taking the curl of the balance equation. The nonlinearity 
of ( |2.1b| ) vanishes in the projection procedure when using only one velocity mode. Now, 
the solution procedure is obvious: Solving the linear system (|3.9| ) for a given convective 
amplitude and inserting the result into ( p.lO|) yields the Rayleigh number of that SOC state 



with the convective amplitude ifmax- Thus, the bifurcation diagram -R(wmax) or tymax(-R) 
may be calculated. 
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Finding the TW solution is slightly more difficult since the modes are time dependent: 



■l)ujt 



The form of the time dependence is determined by the fact that the TWs are stationary in 
a frame comoving with v = uj/k. Therefore, the angular frequency of the nth lateral Fourier 
mode (lower index of the mode amplitudes) is nuj with uj being the basic frequency of the 
TW. With the complex vector of TW mode amplitudes 

T 



X 



TW 



n2(m+l) ;^2m+l 'T^m+l a2' 



- 2n 



2n+l 1 S2n ) S2n+1 



the system of the projected balance equations reduces once more to a linear algebraic system 



) ■ Xtw — ^TwyWn 



(3.11) 



where, however, the matrix ■Mj.^y of mode coupling coefficients is complex. Another complex 
equation (or two real ones) generated by the momentum balance relates the two real groups 
UJ and Wmax with the Rayleigh number R and the separation ratio ip 

( \ ( \ ( \ 



aiiXxw a-i2XTw 



V 



a2iXTW a22XTW 



1/i? 



biX 



TW 



v 



boX 



2-^TW 



(3.12) 



with ajfc, bj being vectors of projection coefficients. 

A possibility to solve the system is to solve (|3.11|) for a given combination (c<j,Wmax) 
and to use the resulting Xtw(i^, "W^max) for solving (|3.12| ) with resepect to {1/R,ip). This 
means that R and if) are uniquely determined for a given combination [uP^w"^^-^. The 
relevant parameters are here the squares u;^ and w^ax since left- and right-traveling waves 
are symmetry degenerated and the convective amplitude Wmax was chosen to be positive by 

(EH)- 

This result is illustrated in Fig. ^: The TW states fill the cu^-w^^x^plane. The lines 
are connecting states along a TW bifurcation branch for a particular separation ratio. The 
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control parameter R varies along a particular line between R^sc at the Hopf bifurcation 
(w^max = 0) and R* at the SOC-TW-transition (cu^ = 0) in a way that can be non- 
monotonous. The value of uj"^ for w^^x = is uj^, i. e., the square of the Hopf frequency 
uh- The dotted line represents the position of the saddle nodes in the diagrams lj{R) or 
u^max(-R)- It vanishes in the vicinity of ip = —0.01 since for weaker Soret couplings only an 
unstable lower TW bifurcation branch exists without any saddle node. These topological 
features will be discussed later in this article. 



D. Comparison with finite difference numerical results 

In Fig. ^ we give a quantitative comparison of bifurcation diagrams computed by the 
above presented numerical Galerkin method and those obtained by a finite difference MAC 
scheme with a full representation of the fields. In our Galerkin scheme, we used one velocity 
mode and up to 19 temperature and ( field modes in each direction. So, we computed 761 
real mode amplitudes in the case of TW convection. 

For both SOC and TW states the pairs of diagrams coincide. The most evident dis- 
crepancy can be observed in the end point r* of traveling wave convection in Fig. 0(d): 
The MAC results with a discretization of dx = dz = 0.05 (solid circles) predict r* ^ 1.65 
|16| , the Galerkin scheme (dashed line) r* ~ 1.495. In order to elucidate this deviation we 



performed a finite difference calculation with dx = dz = 0.025 (open lozenges) and observed 
r* ~ 1.45, i. e., close to the value of the Galerkin scheme which used modes up to a wave 
number 197r in both directions. The variation of r* with the spatial resolution of the MAC 
scheme is caused by concentration boundary layers at the plates: Galerkin method and 
the finer resolving MAC scheme show the same thickness which is smaller than that pre- 
dicted by the worse resolving method dx = dz = 0.05. Now, the SOC-TW-transition may 
be interpreted as a boundary layer instability that occurs if the SOC boundary layer 
thickness oc (L/wraaxY^^ exceeds a certain value when reducing r or reducing the convective 
amplitude. This critical thickness is reached for higher amplitudes when using a numerical 
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method that produces larger boundary layers. In so far, r* ~ 1.45 is a more adequate value 
for L = 0.01, a = 10, and = -0.25 than r* ^ 1.65 



The deviations in the bifurcation diagrams for the variance M of the concentration in 
Fig. |^(b,e) are mainly due to the shift in the frequency bifurcation branches. This means 
that the dependence of the concentration distribution on the frequency is reproduced well. 

The convective heat transport A^ — 1 [Fig. §(c,f)] in our Galerkin approximation is carried 
by only one velocity mode so that the actual values may be expected to be larger the higher 
the forcing, i .e., the Rayleigh number r. This typical behaviour can be observed also for 
the pure fluid = 0. There, a discrepancy between full MAC results and our Galerkin 
approximation of about 6% at r = 2 is seen. 

IV. RESULTS 

A. Soret coupling and bifurcation topology 

The dependence of the bifurcation topology on the strength of the Soret coupling, say in 
the range —0.25 < ip < 0, was not yet discussed in the literature [n,16,2^ in detail: On the 



one hand, the directly integrating numerical methods |T^JTB| require large amounts of CPU 



time due to critical slowing down near the saddle node positions t^qq and r^-^. On the 
other hand, the Soret effect was only implemented incompletely in a theoretical approach 
28| . The same applies to the transition point r* of TW to stationary convection. The 



Galerkin method presented in the last section does not suffer from these disadvantages and 
it yields also the whole unstable branches directly without any numerical tricks. Since it 
computes only flxed points the intrinsic time scale of the system does not enter the problem. 
However, an additional stability analysis of the computed states is necessary. 
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1. Strong negative ■0 

The investigation of the variation of the bifurcation scenario with the separation ratio 



ijj has shown for strong negative couphngs a very interesting feature |17|: the development 
of a lower stable TW branch out of the bump seen in Fig. |l| in the unstable branch for 
if) = —0.25. For ip < —0.4 two stable, convective TW solutions exist opening up the 
possibility of the coexistence of two different traveling states in one container. Furthermore, 
all these TWs were found to display universal scaling properties [[T^: (i) The mixing M, 
i. e., the concentration variance varies linearly with the TW frequency, (ii) The latter itself 
is uniquely fixed by the "distance", r{wjaax,4') ~ "'^('"^maxj = 0), of the TW state's location 
Wmax) on the TW bifurcation branch from the pure fluid {ip = 0) convection coordinates 
in the r-w^^x bifurcation diagram. Thus, bifurcation properties, spatial structure of the 
concentration distribution, and TW dynamics show for strong Soret coupling a peculiar 
scaling behaviour. For a discussion of the characteristic changes in the bifurcation topology 



at strong Soret coupling see Ref. |17 . 



2. Weak negative ip 

On the side of weak negative Soret couplings, i. e., ip ^ 0, the motion of the saddle nodes 
Tgoc and r^^v ^"^^ of SOC-TW merging point r* in the r-tp plane was not elucidated 
except for the vanishing of stable TW convection for ip > —0.01 in mixtures with ethanol- 
water parameters L = 0.01 and a = 10 [l^. To fill this gap we discuss in Figs. the 



bifurcation properties of TWs and SOCs in the range —0.25 < < 0. In Fig. ^ we show 
TW and SOC bifurcation diagrams of if^ax ''^^ several ip as indicated. In the case of 
the strongest coupling ip = —0.05 [Fig. ||(a)] the same situation as for ip = —0.25 (Fig. |l]) 
is observed except for the fact that here r* < rose so that a SOC state is observed when 
heating above threshold. The stationary bifurcation threshold, rgtat, is negative as it is the 
case for all ip < ji^ = fo^ = 0-01 pO|JT9[| : Thus, the SOC branch is disconnected 
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with the ground state w^ax = at the positive r-axis. The shape of the SOC branch does 
not change when reducing the strength of the Soret effect because the tricritical separation 
ratio for SOCs, ^/'soc' scales with —Li^ ||2^ and is effectively for small Lewis numbers. The 



most evident effect is the motion of the SOC-TW-transition point r* along the SOC branch 
towards the SOC saddle node. For ip = -0.02 [Fig. |(b)], the TW and SOC saddle nodes 
coincide; for ip = —0.01 [Fig. |^(c)] the TW branch merges with the SOC branch at the 
SOC saddle node. The transition from SOC to a pair of symmetry degenerated TWs at r* 
which is a pitch fork bifurcation of the TW frequency is still backwards in the sense that the 
frequency bifurcation is subcritical. Different to = —0.02 we can observe for ip = —0.01 
TWs also for control parameters smaller than those of all SOCs, i. e., the extended state 
with the smallest Rayleigh number is now a TW and no longer a SOC. This is different from 
the behaviour for ip < —0.02 where SOCs exist also at Rayleigh numbers not allowing TWs. 

In Fig. ^(d) with ip = —0.0085 the TW branch merges with the lower, unstable SOC 
branch which then becomes stable beyond its saddle node. But the TW saddle is still found 
at smaller amplitudes and Rayleigh numbers than that of the SOCs. The consequences are 
as follows: The lower TW branch locates unstable TWs which become stable at the saddle. 
However, they change stability once more since stable TWs may not merge with the unstable 
SOC branch. The only possibility is the existence of an additional bifurcation point — here 
to a modulated TW (MTW) — on the upper TW branch between and r*. This scenario 



is investigated in more detail in the next subsection [IV A3 . 

Here, we continue the discussion of the changes in the bifurcation topology of ethanol- 
water mixtures for ip ^ 0. At ip = —0.007 in Fig. ^(e) the TW saddle has vanished and 
the whole TW branch is unstable. The next significant change in the topology occurs at 
the tricritical separation ratio ^/'xw — "~5-10^^ [^. Therefore, at ip = — 4-10~^ we observe 
in Fig. ^(f) a supercritically bifurcating TW branch. It is stable just at the onset and then 
becomes unstable with respect to MTWs. At the codimension-2-point ipcT = -3.526-10"^ 
PP| the TW branch vanishes completely. For separation ratios ip between '^/'tw "^ct 
there is also a necessity of a change in stability along the supercritically bifurcating TW 
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branch ending on the unstable SOC branch. 



3. Modulated TWs 

We have checked the scenario for the appearence of MTWs more exphcitly for a parameter 
combination which is reahzed in ^He-'^He-mixtures rather than in ethanol-water, namely 
L = 0.03, a = 1, and ip = —0.055. In ^He-'^He this scenario occurs in a four times broader 
range of the control parameter (r G [1.0612,1.0614], see Fig. P(a)) than it is realized in 
ethanol-water mixtures (r G [1.01375, 1.01380], see Fig. Kd)). 

Sufficiently above the oscillatory threshold rose = 1-08827 {ujh = 3.8467, = 3.1152) 
stable SOCs (filled triangles) can be realized according to Fig. |](a). When reducing the 
heating rate r we approach the SOC-saddle at about t^q^ ^ 1.06139. Below that, the 
system falls down on the branch of TWs which is oscillatory unstable (open squares) at 
rgoc- The upper stable TW branch (filled circles) turns unstable at tmtw — 1-06125 shortly 
above the saddle r^w- The oscillatory unstable TWs (open squares) get modulated by a 
frequency at least ten times smaller than the basic frequency of the TW itself — compare 
the frequencies in Fig. |](b) and the imaginary parts of the relevant eigenvalues in Fig. ^(c). 
The instability of the TWs on their upper branch, i. e., the transition at tmtw from filled 
circles to open squares, gives birth to MTWs with slowly breathing amplitude. The type of 
the instability of the TWs changes on the upper TW branch at r ^ 1.06141 from oscillatory 
to stationary (open circles) in order to merge appropriately with the stationary unstable 
lower SOC branch (open triangles). 

This scenario has also been discussed in a minimal model by Knobloch and Moore [^ . 



They, however, investigated a system with stress free and, more importantly, permeable 
boundaries without adequately resolving the boundary layers. We give here evidence for 
MTWs in binary mixtures with realistic boundary conditions realized in the experiments. 
Note, however, the experimental investigation requires a control of the temperature differ- 
ence of about 10~^. 
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4- Phase diagram 

Our results are summarized in a phase diagram in Fig. |^ where the ^/'-dependence of r* 
(open lozenges), SOC saddle (open triangles up), and TW saddle (filled circles) are displayed 
together with the linear stability thresholds of the basic state (stationary: open triangles 
down, oscillatory: filled squares). Stable TW states on the upper TW bifurcation branch are 
located in the shaded region of Fig. 0. The inset in the upper right corner covers separation 



ratios between ^p = —0.05 and ijj = —0.005. In this inset we have scaled r — 1 by in 
order to map the SOC saddle node position rgoc approximately onto a constant since its 
^/'-dependence may be fitted very well by 



SOC 



1 + 1.636^^1?/ . (4.1) 



However, in view of the fact that the exponent in ( [l.lj) is not exactly 1/2 it is not surprising 
that r^QQ is not exactly constant in the inset of Fig. 0. The range of existence of the MTWs 
is too small to be visible in Fig. ^ 

The changes of the bifurcation topology induced by increasing the strength of the 
negative Soret coupling beyond values of —0.25 including a detailed phase diagram for 



^jj E [—0.7, —0.2] were discussed in Ref. ||T^: There, two stable, nonlinear TW branches have 



been discovered for ip < —0.4 in mixtures with L = 0(0.01) and a = 0(10). 



B. Concentration distribution and streamfunction 

Fig. 1^ shows a combined bifurcation diagram of convective amplitude Wmax and phase 
velocity v vs. reduced Rayleigh number r. The TW states have been computed with our 
Galerkin method using the one- mode velocity field approximation ( p.l|) . Four states are 
labelled by A to D in agreement with Fig. |I[ 

We want to discuss here the changes in the concentration distribution when moving along 
the TW bifurcation branch and their relation to the structural changes of the streamfunction 
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4>{x,z) = sin kxCi{z) + VZ (4.2) 

for the velocity field in the frame that is comoving with the TW phase velocity v. Note 
that in this frame the velocity field is stationary so that passive particles would move along 
the streamlines of ( |4.2| ). 

The first occurrence of local extrema in z) ( [4.2| ) gives rise to the first appearence of 
areas of closed streamlines. That happens for A; = vr at the velocity ratio x = ^^^if^ = 0.9807. 
Thus, for X < 1, i. e., in all TW states between the Hopf bifurcation threshold rose and state 
A the streamlines are all open. Increasing the convective amplitude Wmax yields decreasing 
phase velocity v and increasing values of x so that areas of closed streamlines appear that 
grow on cost of those occupied by open streamlines. For f = 0, i. e., at the point D only 
closed streamlines are observed. 

Without feedback into buoyancy concentration is a passive scalar transported by means 
of convection and diffusion (for the discussion of the Soret coupling in the bulk equations see 
Sec. [IV C|) . Then one can apply the results of passive scalar theory, e. g., [0,0 . It explains 



that within closed streamlines a weakly diffusing scalar is homogenized. This behaviour is 
elucidated in Fig. ^ where we show lateral and vertical concentration profiles for the states 
A - D labelled in Fig. |^. In the states B, C, and D where x > 1 therefore areas of open 
streamlines exist we see a plateau characteristic of both profiles. However, the state A with 
only open streamlines shows no such plateaus. 

Note that the lateral concentration wave profiles between the Hopf bifurcation and the 
state A are basically harmonic and that their amplitudes increase with increasing convective 
amplitude. This amplitude growth of the harmonic concentration wave occurs as long as 
there are not yet closed streamlines, i. e., up to a limiting point in the vicinity of state A. 
The amplitude of the wave profile of the state B in Fig. ^ is larger than that of state A as 
the lateral profiles are not taken at the center of the closed streamlines but in the center 
of the convection cell. Beyond this maximum amplitude the concentration wave crests are 
"cut off" and a plateau develops with the appearence of closed streamlines. The plateau 
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extension, i .e., the region of constant concentration within the areas of closed streamhnes 
broadens with increasing x while simultaneously the plateau value of the concentration wave 
reduces. Furthermore, a small concentration peak at the leading front of the right traveling 
wave develops. It is caused by the advective injection of the concentration at the plates 
into the areas of closed streamlines. The lateral profile of state A has only a very weak, 
not visible asymmetry with respect to reflection at x = ±1/2 whereas the asymmetry in 
state B is obvious and gets more and more pronounced in C and D. The adevctive injection 
of concentration into the regions of closed streamlines takes — with diffusion being small 
compared to advection — the form of an inwards spiralling concentration jet that can be 



seen in Fig. 3 of Ref. ||T6|. 

This strong relation between streamlines and concentration distribution is also demon- 



strated in Fig. |T0| where we have plotted the concentration field C vs. the streamf unction 
in the frame comoving with the TW state traveling to the right. Just at the Hopf bifur- 
cation threshold the streamfunction in the comoving frame is (f) = vz = Unz/k reflecting 
a vanishing convective velocity field in the laboratory frame. Via the conductive concen- 
tration profile Ccond = —"ipz we get the relation Ccond{<t>) = ~^4> — 0.0700 (dotted lines 
in Fig. |r^) between concentration and streamfunction with varying between —v/2 and 
f /2, i. e., between —1.78 and 1.78. In state A, the conductive concentration distribution is 
slightly deformed and the streamfunction varies between —1.4 and 1.4 which are the values 
at the top and bottom plates. At these boundaries the concentration is in contrast to the 
streamfunction not constant so that a vertical shape of C(0) is observed there. In state 
B there are already small areas of closed streamlines with equilibrated concentration. In 
these regions C(0) is a constant. Moving along the TW bifurcation branch the areas of 
these regions keep on increasing via state C to state D which is a SOC state. In this state 
we observe a remarkable concentration variation around its mean value only in a small 
streamfunction interval around 0, the separatrix in this state. This is a strong boundary 
layer phenomenon caused by the smallness of Lewis number L = 0(0.01) in comparison 
with the convective amplitude Wmax = 0(10). 
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C. Soret coupling and concentration current 

1. Comoving frame of reference 

The impact of the Soret couphng on the concentration current in a TW is best understood 
by studying the current in a frame comoving with the TW's phase velocity v. In this frame, 
the velocity field u(x, z) = {—dz(p, 0, dx4>) is well approximated by the ansatz ( [4 .21 ) according 



to ansatz ( |3.1D and the discussion of Fig. 0. The corresponding streamlines can be seen in 
Fig. |lT](a) for a TW state with v = A/tt, i. e., on = 4. 

In the comoving frame the concentration current J is given by 

J = uC - LV{C -ipT) . (4.3) 

In a relaxed TW state the relation 

= = - V ■ J 

holds so that J is divergence free. It is shown as a vector field plot in Fig. |ll](a) in the left 
column together with the streamlines of the velocity field u for a TW propagating to the 
right. The structure of the current can be explained as follows: In the closed streamlines 



of the right rolls of Fig. [TT] concentration is homogeneized to a level close to that of the 
upper plate. This is due to the vicinity of the right roll to this plate where the higher 
concentration (C > 0) is observed for negative Soret couplings and heating from below. 
The left area of closed streamlines is located next to the lower plate and contains therefore 
the lower concentration [C < 0). The two different signs of concentration — the mean 
concentration is normalized to zero — are the reason for the concentration current J to 
rotate clockwise in both clockwise and counterclockwise rotating rolls. In the vicinity of the 
center open streamline that meanders between the rolls, a line of vanshing concentration 
current exists. It also vanishes near the center of the rolls where the velocity u vanishes. 

At the plates the concentration current is purely lateral since due to the impermeable 
boundary condition J ■ = 0. It is mainly given by —vC being negative (positive) at 
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the top (bottom) plate with C > (C < 0). The contribution from lateral concentration 
gradients along the plates to the concentration current ( |4.3| ) is mutiplied by the Lewis number 
L = 0(0.01) and can therefrore be neglected there. The temperature, and consequently the 
Soret effect, does not contribute at all to J • e^, at the plates since the temperature is fixed 
there. 



2. Approximate Soret induced current 

Next, we discuss the influence of the temperature field on the concentration current ( |4.3|) 
in the bulk. We do this by studying the left column of Fig. |ll|(b). There, Lip'VT is plotted 
by arrows whose lengths are magnified by a factor 30 relative to that of J in Fig. |Tl](a). 
For negative Soret coupling Lip'VT is parallel to the diffusive heat current — VT pointing 
upwards in the system heated from below. The important thing is the existence of this 
current and its mean upwards direction and not the small modulations. We ignore these 
small lateral modulations by replacing Ltp'VT by an adequate mean. For this mean we 
choose the mean of Lijj'VT at the plates, namely 

in order to guarantee the impermeability of the plates in the lateral mean. A global averaging 
of Lip'VT would lead to —LipSz differing from our choice only by a factor = 0(1) and 
violating the impermeability by the same amount. 

The replacement of Lip'VT by the mean (|4.4| ) leads to the modified concentration current 



Jmodificd = uC- LVC - LijNe, (4.5) 



shown in the right column of Fig. |TT](a). Therefore, in the modified concentration balance 
= dtC = -V- J^odified = -V • [uC - LVC] (4.6) 

the tempreature field has disappeared. It occurs only in the boundary condition 
= ■ Jmodified = -LdzC - LipN at z = ±1/2 . 
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This allows the concentration field to be described by 

oo oo 

C{x, z; t) = -iIjNz + ^Y^ [C22™+^(t)e-2™'=" + c.c] ^2 sin(2m + l)nz 

n=0 m=0 

oo oo 

+ 5Z 5Z [^2^+1(^)6"^^^"+^^*^'' + c.c] y2cos2m7rz . 

n=0 m=0 

Since now the actual concentration balance is solved and not the balance for the ^-field 
the above discussed problems with the relevance of temperature modes that are 1/L times 
smaller than the ^-field modes do no longer occur. Consequently, the temperature field can 
now be represented by the simplest approximation 

T{x,z,t) = -z + T^{t)V2sm2nz 

+ [T^{t)e~''''' + c.c] V2cos7r2 (4.7) 



just like in the classical Lorenz model |]53|. The advantage of this procedure is that only for 
the concentration field a many mode Galerkin representation is necessary and the tempera- 
ture can be modelled by those modes which are observed in the fields and not additionally 
by those becoming necessary out of numerical reasons. 

Bifurcation diagrams of the frequency computed by this approximate method (solid lines) 



are displayed in Fig. |T2| and compared with the "exact" results (symbols) as described in 
Sec. pi C| over a wide range of Soret couplings. We have chosen bifurcation diagrams for the 
frequency since it is most sensitive to an insufficient mode truncation. For separation ratios 
ip > —0.15 no differences are observed whereas for stronger couplings {ip = —0.25) at high 
heating rates deviations become visible. The reason for this is the neglect of temperature 
modes higher than those incorporated in ( [4.7| ). At high heating rates they become relevant: 
The maximum of the second lateral temperature mode has reached for r = 1.6 about 6% of 



the size of the first one (3% for the third relative to the first) (see also ||16| , Fig.5b]). In so 
far, the simplest ansatz for the temperature field ( [4.7|) cuts off at least 10% of the spectrum 
leading to comparable errors in the bifurcation diagram, especially in the vicinity of the 
SOC-TW-transition. 
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V. CONCLUSION 

We have given a detailed analysis of the influence of the Soret effect on thermal convection 
in binary liquid mixtures. As a tool we used a many mode Galerkin method combined with 
an approximation in the velocity fleld: It is truncated by a single mode which holds for 
convection in fluids with Prandtl numbers cr > 1. The mean flow in traveling waves (TWs) 
was recognized as unimportant since it is very small and contributes only unsystematically 
to the main TW properties. This truncation of the velocity lead to a very simple and 
efficient solution procedure for the nonlinear TW and SOC fixed points by only solving 
linear equations. This opened up the possibility for a detailed elucidation of the changes in 
the combined bifurcation topology of stationary and traveling states, especially the existence 
range of TWs in the control and fiuid parameter plane. Together with the investigation of 



strong Soret couplings interesting fields for experiments are opened: the bistability 
of slow and fast traveling waves and the occurrence of modulated traveling waves with 
slowly breathing amplitude. The classification of the states is facilitated by a detailed phase 
diagram. As an additional insight it was found that the order paramters TW frequency and 
TW convective velocity determine the control parameters Rayleigh number and separation 
ratio uniquely. 

Eliminating the intrinsic time scale from the computation of the fixed points allowed be- 
sides the determination of bifurcation points and saddle nodes the investigation of unstable 
states. In this regime the transition from weakly to strongly nonlinear TWs is observed and 
may be understood in the framework of the intimate relation between concentration distri- 
bution and the structure of the fiow and the changes in this relation along the bifurcation 
branch. 

As a further result the reason for the failure of earlier Galerkin approximations for the 
convection in binary liquid mixtures was revealed: In order to exactly fulfill the concentration 
boundary condition which is coupled to the temperature field earlier approaches used a 
combination of concentration and temperature field without resolving the temperature field 
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adequately. When using this combined field for small Lewis numbers L ^ 1 a resolution of 
the temperature is required that goes beyond the dominating modes seen in simulations and 
experiments. Although the contribution of higher modes to the temperature field are small 
they are essential in the balance equations formulated with the combined field. A way out of 
this artificial theoretical dilemma is obtained by investigating the concentration current: It 
is mainly infiuenced by the lateral average of the temperature gradient in the system. This 
allows to ignore the Soret effect in the bulk equation of the concentration balance (bulk 
Soret effect) and to truncate the concentration field with an adequate boundary condition 
(boundary Soret effect) directly. Then, the foundations of a more compact description and 



solution for the convection in binary liquid mixtures are laid |34 
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FIGURES 

FIG. 1. Bifurcation diagrams of convection flow intensity (a) and frequency (b) in TW (solid 
lines) and SOC (dashed lines) solutions in a binary mixture with ip = —0.25, L = 0.01 and a = 10. 
States A-D are identified for later discussion related to Figs. |-10. 



FIG. 2. Quality of the one-mode velocity field approximation (3J-) for TWs with L = 0.01, 
o" = 10 and ip = —0.25 (left column) and ip = —0.6 (right column) as a function of growing flow 
amplitude it'max along the TW bifurcation branches. The top row displays the error in the vertical 
velocity field according to definition ( p.2| ) and the bottom row the maximal amplitude of the lateral 
mean flow. 

FIG. 3. TW states in the w'^^^-lo'^ -plane for L = 0.01 and a = 10. The separation ratio is 
constant on each solid line. It varies logarithmically form line to line and has the value shown at 
the left ordinate. The dotted line gives the positions of the TW saddle nodes. 

FIG. 4. Comparison of bifurcation diagrams of frequencies uj (a,d), variance M of the concen- 
tration field (b,e), and Nusselt number (c,f) calculated by our Galerkin method and a finite 
difference MAC scheme. Results of the Galerkin method are shown by solid (SOC) and dashed 
(TW) lines. MAC results with a spatial resolution of dx = dz = 0.05 [16] are shown by filled circles 
(TW), open squares (SOC), and open triangles (phase fixed, unstable SOCs). More accurate MAC 
results {dx = dz = 0.025) are displayed as open lozenges. The dotted lines show results of a pure 
fluid, ■0 = 0: thick dots - Galerkin, thin dots - MAC. 

FIG. 5. Bifurcation diagrams of li'max versus reduced Rayleigh number r at weak Soret couplings 
and L = 0.01, a = 10. Open lozenges represent SOCs, filled lozenges TWs. Results were computed 
with our Galerkin method. 
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FIG. 6. Bifurcation diagrams of convective amplitude Wma.x (a) and frequency uj (b) for fluids 
with L = 0.03, a = 1, and ip = —0.055. Graph (c) shows eigenvalues determining the stability of 
the TWs: open (filled) symbols correspond to unstable, Re7Tw > 0, (stable, Re7Tw < 0) states 
(stationary: Im7Tw = 0, oscillatory: Im7Tw 7^ 0). The states corresponding to the small filled 
circles are not shown in (a) and (b). 

FIG. 7. Phase diagram of the ^'"dependence in mixtures with L = 0.01 and o" = 10 in a double 
logarithmical plot r — 1 vs. tp. SOC properties are shown by dotted lines (saddle node r^Q^: open 
triangles up, bifurcation rgtat^ open triangles down). Solid lines correspond to points in the TW 
bifurcation diagrams (Hopf bifurcation rose: filled squares, SOC-TW-transition r*: open lozenges, 
saddle node r^^: filled circles). 

FIG. 8. Bifurcation diag ram of convGctiv6 amplitucic tt^max 

and phase velocity v vs. reduced 
Rayleigh number r for TWs in a mixture with L = 0.01, a = 10, tp = —0.25, and k = n. The 
stable (unstable) branches are shown by solid (dotted for f^max and dashed for v) lines. Letters A 
- D identify states discussed in the text. 

FIG. 9. Lateral and vertical concentration profiles of the four states identified in Figs. [l| and ^. 

FIG. 10. Concentration C vs. streamfunction cf) in the comoving frame for the states identified 
in Figs, m and ^. The dotted lines referring to the conductive state are explained in the text. 
Since C and (f) was evaluated on a grid some of the fine structure seen in the plots reflects the grid 
discretization. 
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FIG. 11. Concentration currents in a TW propagating to the right in the comoving frame 
(L = 0.01, a = W, ip = —0.25, and to = 4) with complete Soret coupling (left column) and 
averaged, boundary Soret forcing (right column). For details see text. The top row (a) shows the 
total concentration current J ( [4. 3D in the left column and the modified current Jmodified (^^) in 
the right column each together with the streamlines of u. The scaling factor for the lengths of the 
arrows is 0.5 relative to the units of lateral and vertical axes. Row (b) shows the influence of the 
Soret coupling more explicitly. Left column: Lip'VT (complete coupling); right column: —LipNez 
(averaged coupling). In row (b), the arrow lengths had to be magnified by a factor 30 relative to 
those of the top row in order to make them visible. 

FIG. 12. Bifurcation diagrams of TW frequencies. Symbols refer to correct Soret coupling and 
full mode field representation. Lines were obtained by the laterally averaged Soret effect ( |4.5D and 
the reduced temperature representation ( [4. 7]) . The unstable states for the complete coupling have 
been dropped for reasons of clarity. Parameters are L = 0.01 and a = 10. 
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Fig. 10: Hollinger and Lucke, Phys. Rev. E 
Influence of the Soret effect on convection of binary fluids 
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